Electric Dipolar Susceptibility of the Anderson-Holstein Model 
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The temperature dependence of electric dipolar susceptibility \p is discussed on the basis of the 
Anderson-Holstein model with the use of a numerical renormalization group (NRG) technique. Note 
that \p is related with phonon Green's function D. In order to obtain correct temperature depen- 
dence of xp at l° w temperatures, we propose a method to evaluate \p through the Dyson equation 
from charge susceptibility Xc calculated by the NRG, in contrast to the direct NRG calculation of 
D. We find that the irreducible charge susceptibility estimated from \ c agree with the perturbation 
calculation, suggesting that our method works well. 

PACS numbers: 75.20.Hr, 71.38.-k, 75.40.Cx 

Keywords: Electric dipolar susceptibility, charge susceptibility, numerical renormalization group method 



I 

O 

o 



>: 

ov 
ov 
in 

o 



In the research field of condensed matter physics, ex- 
otic magnetism in cage structure materials such as filled 
skutterudites has attracted much attention due to the 
interests on new electronic properties caused by oscilla- 
tion of a guest atom in a cage composed of relatively 
light atoms Such oscillation with large amplitude is 
frequently called rattling and it is considered to play cru- 
cial roles for the formation of magnetically-robust heavy 
electron state in SmOs4Sbi2 [2J. This peculiar heavy- 
electron state has been theoretically investigated from 
various aspects by several groups [3T-[l7l]. 

Recently, the Kondo effect of a vibrating magnetic ion 
in a cage has been theoretically discussed on the basis 
of a two-channel conduction electron system hybridized 
with a vibrating magnetic ion [18-20]. Note here that a 
vibrating ion inevitably induces electric dipole moment. 
Then, it has been found that magnetic and non-magnetic 
Kondo effects alternatively occur due to the screening of 
spin moment and electric dipole moment of vibrating ion 
|2l| . In particular, electric dipolar two-channel Kondo 
effect has been found to occur for weak Coulomb inter- 
action. Then, it has been proposed that magnetically 
robust heavy-electron state appears near the fixed point 
of electric dipolar two-channel Kondo effect. 

In this paper, in order to promote our understanding 
on the Kondo effect concerning electric dipole moment 
P, we analyze the temperature dependence of electric 
dipolar susceptibility \p on the basis of the Anderson- 
Holstein Hamiltonian with the use of a numerical renor- 
malization group method. For the reproduction of cor- 
rect temperature dependence of \p at l° w temperatures, 
we propose a method to evaluate xp through the Dyson 
equation from charge susceptibility Xc- This method is 
found to provide correct results in the temperature region 
lower than the Kondo temperature, in sharp contrast to 
the numerical evaluation of the phonon Green's function 
which is directly related to xp- 

Now we explain the model Hamiltonian. We consider 
a conduction electron system in which an impurity ion 



is embedded. On the impurity site, localized electrons 
are coupled with ion vibration. The situation is well de- 
scribed by the Anderson-Holstein model, given by [22[ 
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where is the dispersion of conduction electron, Cka 
is an annihilation operator of conduction electron with 
momentum k and spin er, d a denotes an annihilation op- 
erator of the localized electron with spin a, and V is 
the hybridization between conduction and localized elec- 
trons. The local term H\ ac is given by 



H\ oc = tip + gxp + p 2 /2 + uj 2 x 2 /2, 



(2) 
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where p is a chemical potential, p — ^ a n„ with n a — 
d\d a , g denotes the coupling between electron density 
and ion vibration, x is normal coordinate of the vibrating 
ion, p indicates the corresponding canonical momentum, 
and lu is the vibration energy. Note that the reduced 
mass of the vibrating ion is set as unity. 

In the present model, we ignore the Coulomb interac- 
tion term Un^n^. Of course, the effect of U is quite im- 
portant, but our main purpose here is to understand the 
temperature dependence of electric dipolar susceptibility 
Xp- Since xp should be affected by U, it is necessary to 
grasp in the first place the temperature dependence of 
Xp without the effect of U. Thus, in the present work, 
the Coulomb interaction is ignored. The effect of U will 
be discussed elsewhere in future. 

For actual calculations, it is convenient to introduce 
phonon operator b and b^ through the relation of x — 
(b + tf)/y/2cj. Then, the local term is rewritten as 

Hioc =pp + V^ojp{b + tf) + uj{b ] b + 1/2), (3) 

where a is the non-dimensional electron-phonon cou- 
pling, defined by a — g 2 /(2uj 3 ). Concerning the average 
electron number, throughout this paper, we consider the 
half-filling case, at which p is given by p — 2au. 

The conduction electron model hybridized with local 
impurity is precisely analyzed with the use of a nu- 
merical renormalization group (NRG) technique [23l.l2ij]. 
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FIG. 1: (Color online) (a) Entropy 5j mp and specific heat Ci mp 
vs. temperature T for a = 2 and lj = 0.2. (b) Charge sus- 
ceptibility T\c and spin susceptibility T\s vs. temperature 
T for the same parameters as (a). 



The logarithmic discretization of the momentum space 
is characterized by a parameter A and we keep M low- 
energy states for each renormalization step. Throughout 
this paper, Here we set A=2.5 and M=5000. The energy 
unit is a half-bandwidth of the conduction electron. In 
this unit, we fix V as V — 0.25. As for the calculation of 
phonon part, the number of phonon basis is 50. 

In order to clarify electronic properties of H, first we 
evaluate entropy Si mp and specific heat Ci mp , which are 
given by S ilap =-dF/dT and C imp =-Td 2 F/dT 2 , respec- 
tively, where F is the free energy of local electron and T is 
a temperature, defined by T = A - ^ -1 ^ 2 with the renor- 
malization step number N. We also evaluate charge and 
spin susceptibilities, Xc and Xsi which are, respectively, 
given by 

1 ^ e" Bi / T - e~ E i/ T 
Xc = 7 E 5 5 \(*\P - WW > W 
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and 



Xs= Z^ P~ E, l(i|Sz|j)l 



(•5) 



where £^ is the eigen-energy of H, \i) is the corresponding 
eigen-state, Z is the partition function given by Z = 
£ ie ^/ T , (p) = Y.S\p\^- Ei,T /Z, and Sz =nt - n ; . 

Now we move on to the NRG results. First we show 
typical results for the Kondo phenomenon due to the 
electron-phonon coupling. In Fig. HJa), we depict Simp 
and Ci mp as functions of T for a = 2 and w = 0.2. In 
the high-T region as T > 10~ 2 , 5i mp shows a large value 
due to phonon excited states. For 10~ 7 < T < 10 -2 , 
Simp indicates a plateau of In 2 and C- lmp becomes almost 
zero in the corresponding temperature region. Note that 



we do not understand what degree of freedom is relevant 
to In 2 only from the results of Si mp and Cj mp . This 
point will be discussed later. Around at T = 1.8 x 10~ 9 , 
we observe the release of entropy In 2 and the specific 
heat Ci m p forms a clear peak, which defines the Kondo 
temperature Tk- Thus, we obtain Ik = 1.8 x 10~ 9 for 
the present parameters. Finally, for T < 10~ 9 , both Si mp 
and Cmp vanish, indicating that the system is in the local 
Fermi-liquid state. 

In order to clarify the degree of freedom relevant to 
the present Kondo phenomenon, we evaluate charge and 
spin susceptibilities. The results are shown in Fig. Hlb). 
We immediately notice that the spin susceptibility \s is 
rapidly suppressed even at high temperatures as T ~ 
10 _1 , when we decrease the temperature. On the other 
hand, the charge susceptibility Xc is rather increased to 
the enhanced value for T ~ 10~ 2 . In the temperature 
range of 10 -7 < T < 10~ 2 in which In 2 plateau appears, 
we find the large value of Xc- Around at T = Xk, we 
observe that Xc is gradually suppressed. Namely, the 
Kondo behavior clearly appears in the temperature de- 
pendence of Xc- We conclude that in the present model, 
the charge Kondo effect occurs, since degenerate vacant 
and double occupied states play roles of pseudo spins. 

As for the understanding of the charge Kondo effect, 
it seems to be enough to evaluate Xc- However, in order 
to visualize the situation of the charge Kondo effect, it 
is useful to recall that the vibrating ion induces electric 
dipole moment P, given by P = zex, where z denotes 
the valence of the guest ion, e indicates electron charge, 
and x is ion displacement. Namely, the electric dipole 
susceptibility xp is also related to the Kondo effect for 
the vibrating ion problem [2l| . As easily understood from 
the definition of P, we obtain xp as 
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where D denotes the zero-energy component of the 
phonon Green's function, given by 



1 ^ er E '/ T 

D ^=Z^ E—E~ 

ij 



-E 3 /T 



\(i\u-(u)\j)\ 2 .(7) 



Here u is given by u = b + and we add a subscription 
"NRG" to show explicitly that this quantity is evaluated 
by the NRG method. 

Of course, we can perform the NRG calculation for 
-DnrGj but we consider an alternative way to evaluate 
Xp without using P^nrg- For the purpose, we exploit the 
Dyson equation, which relates Xc and D. In general, the 
Dyson equations are diagrammatically shown in Fig. [2] 
and the first one (a) is expressed by 



D(iv n ) = D (iv n ) +g 2 D (iv n ) U(iu n ), 



(8) 



where v n = 2?:Tn is the boson Matsubara frequency with 
an integer n, D is the dressed phonon Green's function, 
g is the electron-phonon coupling constant, n denotes 
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FIG. 2: Dyson equations for phonon Green's function D 
(thick wavy line) and polarization function IT (bubble with 
hatched square). Thin wavy line and bubble with hatched 
triangle denote non-interacting phonon Green's function Do 
and irreducible polarization function EE, respectively. 



polarization function, and Do is a non-interacting phonon 
Green's function, given by Do(iv n ) = — 2w/[(iv n ) 2 — cj 2 ]. 
At the static limit of v n = 0, by noting that 11(0) = Xc, 
-Do(O) = 2/lo, and g — •Jau), we obtain a relation between 
D and Xc as 

D = 2/uj + 2a Xc - (9) 

In principle, D of eq. ((9]) is equal to Dnrg of eq. flJJ). 

In Fig. H we depict TD NRG and TD vs. T for w = 0.2 
for the comparison of Dnrg and D. We show a couple 
of results for a = 1 and 2. In the high-T region, we 
find that TDnrg and TD agree well with each other. 
To be honest, we observe a little quantitative difference 
between TDnrg and TD, although it is difficult to no- 
tice it in the graph of the logarithmic scale. This type of 
small deviation can be overcome by the elevation of the 
numerical accuracy. For instance, if we keep more num- 
bers of states in each renormalization step, the accuracy 
is expected to be improved. 

However, in the low-temperature region of T < Tk, 
where T K = 1.8 x 10~ 9 for a = 2 and T K = 6.6 x 10~ 4 for 
a = 1, we observe serious discrepancy between TDnrg 
and TD. From the analysis of numerical data, we find 
that TD NRG oc T 2 and TD oc T for T < T K . When we 
fix a temperature and increase a, both TDnrg and TD 
are monotonic increasing functions. The slope of TDnrg 
or TD seems to be independent to a. 

This difference seems to be so serious that it is difficult 
to determine which is the correct behavior. However, 
here we recall the local Fermi- liquid theory for T < Tk- 
Namely, we obtain the narrow quasi-particle band at the 
Fermi level in the electron density of states, leading to 
constant density of states po- Since Xc can be related to 
Po at low temperatures, it is natural to consider that D 
should be constant at the same temperature region. If 
we further imagine that phonon frequency U) is decreased 
due to the Kondo screening, the dressed phonon Green's 
function D = 2/ui should be larger than Do(0) = 2/lj. 
We deduce that u/ui is a monotonic decreasing function 
of a. In short, it seems to be natural to consider that 
TD ex T in the charge Kondo effect. 

If so, here we have a naive question concerning the 
reason why we cannot evaluate correctly D by the NRG 
method. In eq. (O, Dnrg can be calculated in the NRG 




FIG. 3: (Color online) D and Dnrg vs. temperature T for 
uj — 0.2. The solid (dashed) curve indicates TDnrg for a = 1 
(2), while the cross (plus) indicates TD for a — 1 (2). 

method, but it includes implicit problems in the calcu- 
lation accuracy. In the NRG method, in order to deal 
with the states which increase rapidly in each renormal- 
ization step, only M low-energy states are kept and the 
rests are simply discarded. In the renormalization pro- 
cess, the electron-phonon excited states which are indis- 
pensable for the description of phonon excitation will be 
lost. This type of problem has been already pointed out 
in the analysis of the two-channel model [21(. On the 
other hand, since the local charge susceptibility Xc is the 
physical quantity with large contribution from the elec- 
tron states near the Fermi level, we expect that Xc is 
obtained in good accuracy with the use of NRG method. 

Let us check the NRG calculation from a different view- 
point by the comparison with the perturbation calcula- 
tion. Here we define the irreducible polarization function 
IT(w n ), which is related with n and D through the Dyson 
equation (b) in Fig. [5] as 

U{iu n ) = U(iu n ) + g 2 fl(iv n ) 2 D(iv n ). (10) 
At v n = 0, by noting eq. © and Xc — n(0), we obtain 

Xc 1 = Xc 1 + 2aw, (11) 
where irreducible charge susceptibility is defined as 

Xc=n(o). 

In Fig. |4j we show Xc vs. a obtained from eq. (fTT]) with 
the use of the NRG result of Xc for ui = 0.2. Note that 
Xc at low enough temperature is estimated at T = 1.1 x 
10~ 10 . For each value of a , we have checked that Xc is 
unchanged even at lower T. For small a, it is shown that 
the NRG results agree quite well with the perturbation 
result in the order of a. For large a, D is enhanced and 
Xc becomes large, indicating that Xc approaches the value 
of 1/t/pb, where U p h denotes the magnitude of attractive 
interaction, given by D p h = 2au> in the present case. The 
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suppressed as J7 p h = 2aw — U. For U > 2au, the charge 
Kondo effect disappears and the spin Kondo effect occurs 
instead [Hj]. We think that it is interesting to consider 
the behavior of \p f° r U > 2au>. 

In summary, we have discussed the temperature de- 
pendence of electric dipolar susceptibility \p 01 the 
Anderson-Holstein model with the use of the NRG tech- 
nique. We consider a direct method to evaluate phonon 
Green's function -Dnrg and another indirect way to cal- 
culate D from Xc through the Dyson equation. After 
careful investigations, we have concluded that D provides 
correct temperature dependence in contrast to Dnrg f° r 
T smaller than the Kondo temperature. The effect of the 
Coulomb interaction on xp is one °f future problems. 



FIG. 4: (Color online) Irreducible charge susceptibility Xc 
(solid circles) vs. a for T = 1.1 X KT 10 and u = 0.2. The 
solid line indicates the perturbation calculation result of \ c 
up to the first order of a and the dashed curve means the 
expected behavior of Xc ~ l/(2acj) for large a, respectively. 

NRG results for Xc reproduce the behavior of l/(7 p h for 
large a. Thus, it is confirmed that Xc is constant at low 
temperature from eq. (fTTj) . 

Finally, let us provide a comment on the effect of 
Coulomb interaction, which has been completely ignored 
in this paper. When we include U, we imagine that C/ p h is 
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